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Transition energies and oscillator strengths of excitons in dependence on magnetic field are in- 
vestigated in type I and II semiconductor nanorings. A slight deviation from circular (concentric) 
shape of the type II nanoring gives a better observability of the Aharonov-Bohm oscillations since 
the ground state is always optically active. Kinetic equations for the exciton occupation are solved 
with acoustic phonon scattering as the major relaxation process, and absorption and luminescence 
spectra are calculated showing deviations from equilibrium. The presence of a non-radiative exciton 
decay leads to a quenching of the integrated photoluminescence with magnetic field. 
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I. INTRODUCTION 

The Aharonov-Bohm effect (ABE) rests upon the ac- 
tion of the vector potential on quantum mechanical par- 
ticles. The idea behind is quite simple: A charge par- 
ticle orbiting around a region with non-zero magnetic 
flux = / B • dS aquires an energy which is a peri- 
odic function of the magnetic flux $b, the period being 
given by the magnetic flux quantum h/e. Shortly after 
its theoretical predictioufi the ABE has been observed 
experimentally.^^ Doped semiconductor nanorings are 
among actual realizations^ where the ABE could be de- 
tected with high precision. 

The exciton Aharonov-Bohm effect is a more recent 
invention: The neutral exciton consisting of electron and 
hole is predicted to show a similar oscillating behavior. 
Originally, a related many-body system has been studied 
by Wendler and coworkers^ who have considered the 
Coulomb effect when placing two electrons into a quan- 
tum ring - instead of electron and hole as in the present 
exciton case. They form a kind of Wigner molecule with 
interesting properties like crystallization. The approxi- 
mations used originally allowed for an analytic solution 
of the ground state and the related persistent current. 
Later on, a full calculation was carried out, but always 
using a rigid lateral confinement in the ringfi 

For the exciton Aharonov-Bohm effect (X-ABE), the 
limiting case of a nanoring with zero width has been stud- 
ied first as well4i^ The simplicity of the model allowed to 
investigate many aspects of X-ABE like the Berry phase 
and persistent current^ absorption^ii or deviation from 
circular symmetry due to the presence of impurities.— 
One basic result was that the X-ABE can only be ob- 
served as long as the ring diameter is comparable or below 
the exciton Bohr radius a^.— Since the ring size cannot 
be made arbitrarily small, a better idea is to increase as, 
i.e. to weaken the exciton. This can be achieved by sep- 
arating electron and hole using a static electric fieldp^ii^ 
or going to type II material combinations where electron 
and hole are confined in different regions i-'^^d^'^^d^ 

In the original ABE idea, the particle path was as- 



sumed to lie in a region with zero magnetic field, which 
means to concentrate the magnetic fiux into the middle 
of the ring. For a semiconductor nanoring of nanome- 
ter size, this is technically very demanding - if not im- 
possible using state-of-the-art techniques. Therefore, for 
X-ABE a situation with the magnetic field penetrating 
homogeneously the nanoring is more realistic. Then, the 
energy in dependence on magnetic field acquires a non- 
oscillating part roughly proportional to as well which 
can complicate the extraction of X-ABE oscillations enor- 
mously. 

First calculations for excitons in a circular nanoring 
with finite width and homogeneous penetrating magnetic 
fieldiSiis^ could not find the X-ABE for the ground 
state. Theoretical progress was made by calculations on 
two-dimensional annular lattices 3i X-ABE oscillations 
were found but the energy shift was not included. Our 
recent calculations for type I and II nanorings with finite 
width in Ref . [l^ have clearly shown that the ground state 
energy has an oscillatory component, which is related to 
the exciton persistent current. For extracting its ampli- 
tude, we have proposed to calculate the second derivative 
of the energy with respect to B. On the search for op- 
timal nanoring parameters, we have discussed in detail 
different material combinations for type II nanorings. 

Compared to Ref. [l^, the present work deals with 
non-circular geometry and gives a detailed study of ex- 
citon kinetics in order to calculate exciton occupation 
which is not necessarily in equilibrium with the lattice 
temperature. We restrict ourselves to incoherent and 
steady-state excitation following Ref. [2^. A first attempt 
to calculate the exciton photoluminescence (PL) in zero 
width nanorings with a single impurity (and also quan- 
tum dots with respect to X-ABE) has been presented 
recently in Ref. [l^. As a simplification, however, an equi- 
librium Maxwell-Boltzmann distribution for the excitons 
has been assumed, and only the integrated PL intensity 
was shown. We have started with a more realistic exciton 
kinetics, still in the zero width model for circular nanor- 
ings in Ref. HI- Non-radiative decay channels have been 
identifled to be responsible for the decrease of the inte- 
grated PL signal with increasing magnetic field, which is 
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FIG. 1: Schematic drawing of the investigated non-circular 
(eccentric) nanoring (shaded) embedded in a quantum well 

named PL quenching. 

The experimental situation for the X-ABE is far from 
clear: An ensemble of InP/GaAs type II quantum dots 
has been studied in Ref . [2^ and a theoretical explanation 
based on Ref. [1^ indicated some X-ABE oscillations in a 
single dot. Later on, the non-circular shape of the quan- 
tum dot {Di and D2 symmetry [see Section II]) has been 
taken into account.— Furthermore, in a recent single dot 
experiment^ on InP/GaAs quantum dots (grown under 
different conditions) no oscillations have been observed. 
This is consistent with our recent calculation for embed- 
ded InP / GaAs quantum dots in Ref. [H. In a nanoring, 
only the ABE for negatively charged excitons - trions - 
has been observedi^ The X-ABE in nanorings still waits 



where upper (lower) sign refer to electron (hole), g* are 
effective g factors, /is = eh/2ma is the Bohr magne- 
ton, and (T^ the Pauli spin matrix (along z). The vector 
potential is used in symmetric gauge, A(r) = x r. 
The Coulomb potential between electron and hole (last 
term) is screened by the static dielectric constant 65 of 
the semiconductor material. 

Angular momentum and spin parts of the wave func- 
tion can be separated for the heavy-hole exciton, form- 
ing a quadruplet with the z-projection of the total an- 
gular momentum M = ±1 (optically active or bright 
states, with circular polarization cr^) and M = ±2 (dark 
states). Combining electron and hole g factors into ef- 
fective exciton g factors^Sii^ as ±1 ~ ^idh + ffe) ^^'^ 



for its unambigous experimental verification. 

The paper is organized as follows. The exciton Hamil- 
tonian is given in Section II in radial and angular vari- 
ables, and its matrix elements are analyzed according to 
the ring symmetry using arguments of group theory. Ex- 
amples for energies, oscillator strengths, and oscillation 
amplitudes of the X-ABE are shown in Section III, for 
both type I and type II structures. Kinetic equations for 
the exciton occupation including acoustic phonon scat- 
tering and several decay channels are discussed in Section 
IV. Its numerical solution allows to plot absorption and 
photoluminescence spectra for various parameters (Sec- 
tion V). The paper is concluded in Section VI, while more 
technical details are deferred to the Appendix. 



II. EXCITON HAMILTONIAN 

For the heavy-hole exciton in nanostructures, effective 
mass theory with appropriate in-plane (ma) and growth 
direction masses (nia^z) is used often. The subscript 
denotes electron (a = e) and hole quantities (a = h). 
Within this approximation, we will investigate a struc- 
ture which consists of a narrow quantum well with con- 
fining potentials Ua{za) into which a nanoring of general 
symmetry is embedded. The corresponding lateral con- 
finement is given by Va(ra). The three-dimensional vec- 
tor r is decomposed into its two-dimensional in-plane part 
and the z-component in growth direction, = (ra,Za). 
The ring structure is schematically plotted in Fig. [T] In- 
cluding a constant i?-field in z-direction (perpendicular 
to the quantum well), we have the single exciton Hamil- 
tonian 
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9x ±2 ~ ^(.9*1 ~ 9e) afiows to write the spin-dependent 
contribution to the exciton energy as 

^m"' = \9*x,M f^B B , (2) 

which gives rise to the Zeeman splitting linear in B. In 
the following this spin contribution is not written explic- 
itly since its addition to the exciton energies is straight- 
forward. The exciton exchange interaction gives rise to 
an additional fine structure splitting which is neglected 
here in view of the dominant i3-field effects. 

Due to the confinement strength in z-direction, we fac- 
torize from the total wave function the product of sub- 
level wave functions^^, of the lowest electron and hole 
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states in the quantum well. The remaining in-plane 
part of the wave function has to be calculated from 
the z-averaged Hamiltonian.— Introducing polar coor- 
dinates and difference and sum angle [0 = (jj^ — (jih, 
<f> = {(f>e -I- 0h)/2], we arrive at 
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where -I- (— ) sign refers refers to electron (hole). The con- 
finement averaged Coulomb potential Vc{re,rh, (f)) has 
been defined in Ref. [l^. 

The exciton eigenfunction of the state a as solution of 



(4) 



is expanded into the basis of angular momentum cigcn- 
functions, 

^Jre,rh, 0, *) = ^ E ^ILAre, rnV^e''^'' . (5) 
The azimuthal boundary conditions have the usual form 
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aire, Th, (f>e, (f>h +2Tr) , (6) 



which leads using Eq. ([5]) to the following relation be- 
tween / and L: 



L even: / integer; L odd: I half integer . 



(7) 



The expansion functions ui^L,a{re,'rf,.) obey a coupled 
system of Schrodinger equations. The confining potential 
produces the matrix elements 
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where again 4- (— ) sign corresponds to the electron (hole) 
From the reahty of Va{raAa) follows V^'^ = {V^f,)* . 

The confining potential can have an arbitrary symme- 
try (see Fig. [5] for examples) which is classified by the 
two-dimensional point group. There are two types: (i) 
Cn which consists of all rotations about the origin by 
multiples of the angle 27r/n and (ii) I?„ which adds to 
the rotations of C„ reflections with respect to n axes 
passing through the origin. Within C„, the potential 
Vail" a, (f'a) IS invariant under rotations by multiples of the 
angle 27r/n. Changing the integration variable in Eq. ([5]) 
by 0^ = </)(j + 27r/n we obtain 
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FIG. 2: Schematic picture of nanorings with different sym- 
metries. 



TABLE I: The symmetry of the ring geometry leads to re- 
strictions for the matrix elements of the in-plane potential 
Eq. ©. 



Point group of the ring 


Non-zero elements of V^" at 
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k = jn {j integer) and Vj^ = V°fe 


0{2) 


fc = 



Comparison with Eq. ([5]) shows that only elements with 
k = jn (j integer) can be non-zero. Using similar ar- 
guments the reflections in £)„ imply ~ to hold. 
The symmetry properties are summarized in Tab. |T1 

In the absence of any symmetry (Ci) all matrix ele- 
ments in Eq. ([5]) can be non-zero. For increased sym- 
metry {Cn or Dn, n > 1) the matrix decomposes into n 
block matrices, since the matrix elements {lL\Va\l' L') are 
non-zero only \i L — L' = n. In the limiting case of large n 
the point groups C„ and D„ converge to the point group 
of the circle 0(2). Due to the rotational invariance of the 
circle, one degree of freedom can be factorized. This is 
the total angular motion since the corresponding commu- 
tator with the Hamiltonian vanishes, [H, —iUd/d^] = 0. 

The Coulomb potential Vc{re,rh,4') depends only on 
the relative angle, and its matrix elements are diagonal 
in L 

{lL\Vc{r„ru,dp)\l' L') = 6ll'V^ {re,rh) , (10) 
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Putting all terms together, the matrix elements of the 
Hamiltonian are given by 
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The azimuthal kinetic terms can be rewritten as 
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with the exciton reduced mass fix ~ "ie'7?/i/(?Tie + m/j). 
Further, rjc is an effective ring radius for the exciton and 
p a phase shift, defined as 
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Equation (|T2|) resembles the zero-width ring model dis- 
cussed in Ref. HH. For L ~ 0, the last line is responsible 
for the overall increase of energies quadratically in _B, 
while the second line is the source of the X-ABE oscilla- 
tions. We will use Eq. p2p in Section III for extracting 
a reasonable estimate of the B oscillation period. 

For dipole-allowed optical interband transitions, the 
oscillator strength of the exciton state a is given by 
the amplitude of finding electron and hole at the same 
placer^ 



fa = dcv J dr^a{r, r) , 



(14) 



where dcv is the interband dipole matrix element. Intro- 
ducing the single sublevel approximation and the expan- 
sion Eq. (I5|), we find 
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Only the wave function component with L = con- 
tributes to the oscillator strength. States with non-zero 
oscillator strength are called bright states, while dark 
states have = 0. 

As an example of a structure with low symmetry, we 
investigate a ring with Di symmetry. Our simple model 
has two circles as boundaries. The inner one is centered 
at the origin (radius ri), while the outer one (radius 
is displaced by 6, as indicated in Fig. [3l Electron en- 
ergy spectra have been investigated for such a structure 
in Ref. [3^ . calling the ring eccentric. We will use in what 
follows non-circular and eccentric as equivalent names, 
contrasting to concentric and circular. The potential is 
assumed to be constant with a value Vq between these 
two rings, and set to zero outside. In the figures, how- 
ever, zero of energy is set to the confinement gap of the 
quantum well made from the ring material. 

The matrix element Eq. ^ reduces to 
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FIG. 3; Schematic picture of a nanoring with D\ symmetry. 
Within the gray region, the potential is nonzero. Inner and 
outer ring boundaries are circles with radii r\ and r2, the 
centers of which are displaced by h. 



The unit step functions determine the integration bound- 
aries to be ±7r or ±<f>[ra) where 



cos{(l>{ra)) 
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The result can be given analytically, 
Vkira) = V^[5u^oe{ra-ri)e{r2~h~ra) 

+ -^sin{k(j>{ra))9{ra-r2 + b)e{r2+b 



(17) 



(18) 



which properly satisfies the relation V^" = ^-k inherent 
to Di symmetry. 

The circular ring is contained as a special case: Setting 
b = 0, only the first line in Eq. contributes. 



Vkira) = Vo^6k,o0ir2 - ra)0(r, - n) . 



(19) 



Consequently, the Hamiltonian matrix is fully diagonal 
in L, and the wave function expansion Eq. ([5]) reduces 
to the sum over I, while L = is fixed and therefore a 
good quantum number. 



III. X-ABE OF THE EXCITON GROUND 
STATE 

We begin with a discussion of the X-ABE for the exci- 
ton ground state in concentric and eccentric nanorings of 
type I and II as schematically plotted in Fig. |4l The ma- 
terial parameters used are summarized in Tab. [ill The 
effective masses are chosen according to the material in 
which the particle is found predominantly. The values of 
the mass density, the sound velocity, and the deformation 
potentials are taken for the ring material. 

In the calculation, i?-field strengths up to _B = 25 T are 
used which can be easily achieved in experiment. The ra- 
dial coordinates in Eq. (jlip have been discretized on a 
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FIG. 4: Schematic view of in-plane geometry (top) and en- 
ergy profiles for conduction and valence band (bottom) for the 
investigated nanorings of type 1 a) and type II b). Specific 
electron and hole positions are visualized. 



TABLE 11: Confinement potentials Vq , static dielectric con- 
stant es, in-plane effective masses rUa, g factors g^, mass den- 
sity pM, sound velocity s, and deformation potentials for con- 
duction (Dc) and valence band (-Di,). 



GaAs/Alo.23Gao.77As InP/Gao.5ilno.49P 



Vo^ (meV) 


-257 " 


-600 ' 


Vo^ (meV) 


-110 " 


50 ' 




12.5 " 


12.6 ' 


me/mo 


0.067 " 


0.077 ' 


mh/mo 


0.36 " 


0.6 


gl 


0.1 


1.6 


9h 


-1.2 3 


-3.0 


PM (kg m"^') 


5370 ' 


4810 " 


s (m/s) 


5330 ^ 


5230 ■= 


(eV) 


7.0 ' 


6.0 " 


D,, (eV ) 


-3.5 ^ 


-0.6 " 



''Ref. 
'^Ref. 
"^Ref. 
"^Ref. 
/Ref. 
9Ref. [33. 
''Ref.E^ 



grid of 40 points with a grid step of 0.5 nm. The expan- 
sion of the wave function into / and L components has 
been checked for convergency of the results, leading to 
a truncation of |^| < 13 and \L\ < 10. The subsequent 
numerical diagonalization was performed with the im- 
proved Lanczos method^ and checked with the Leapfrog 
mcthod4^ii^ 



A. Type I nanoring 

First, a nanoring of GaAs embedded in the surround- 
ing material Alo.23Gao.77As is considered. This type I 
structure has a confinement of both electron and hole 
within the nanoring. 

The dependence of the lowest three exciton energies 
on magnetic field is plotted in Fig. [5^. In the case of 
circular symmetry (6 = nm, black curves) the states 
can be sorted according to their quantum number L^. 
The ground state has La ~ (solid), the first and the 
second excited ones La = I (dashed) and Lq = — 1 (dot- 
ted), respectively. The energies of the first excited states 
are degenerate at B ~ OT. This degeneracy is lifted at 
B with a relatively small splitting since the effec- 
tive radial distance between electron and hole is small. 
In order to see this dependence more clearly, let us look 
upon Eq. ((T2)) . Due to the strong radial confinement, 
the coordinates Va in Eq. (fT2|) can be replaced by their 
expectation values Ra = {ra)- Their values being rather 
close in the present case {Re = 8.0 nm, R^ = 7.8 nm), we 
realize that L = ±1 gives only a minor difference in the 
energy. 

The case of non-circular symmetry (6=1 nm, red 
(grey) curves in Fig.[S^) is qualitatively not different from 
the circular one since the energetically lowest states for 
different quantum numbers L do not cross. Their mix- 
ing leads to a larger splitting among states and lifts the 
degeneracy at i? = T. More important are the changes 
in the oscillator strength (Fig. [5Jd). The Di symmetry 
of the eccentric nanoring implies that there is one sym- 
metry axis, let us say the x-axis, with reflection operator 
Tx . The wave function transforms in the following way 



T^'^a{Xe,ye,Xh,yh) 



±'^a{Xe,ye,Xh,Vh)- (20) 



With respect to this symmetry operation, all states can 
be grouped into even and odd ones. The odd states have 
zero oscillator strength which follows immediately from 
Eq. (|20p . The doubly degenerate states of the concetric 
nanoring at _B = T can form an even and odd lin- 
ear combination with respect to T^. As the symmetry 
is lowered these combinations get mixed. The energetic 
order of even and odd states can be estimated for the 
lowest ones. The ground state is always even as shown 
in Fig. Oj. The first excited is odd since the contribu- 
tion of the Coulomb interaction is still large as in the 
ground state, while the kinetic energy is lower than for 
the next even state. Further excited states have differ- 
ent Coulomb contributions and that is why their order 
cannot be determined in general. 

All lines in Fig. [5^1' seem to shift upwards quadratically 
in B, which is indeed the dominant part of the diamag- 
netic shift. In order to extract the tiny amplitude of the 
X-ABE oscillations, we have proposed in our previous 
publication Ref. [^^ to calculate the second derivative of 
the exciton energy with respect to the S-ficld. This is 
shown in Fig. [S^. 
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FIG. 5: (Color online) Calculated exciton properties in nanorings in dependence on magnetic field. The left column refers to a 
type I system (GaAs/AlGaAs) with ring radii n — 4 ran, r2 = 12 nm, while for the right column a type II system (InP/GalnP) 
with ri = 8nm and r2 — 16 nm has been considered. Results for a circular ring {b — 0) are shown in black, while red (grey) 
curves are obtained for a non-circular ring (type I: 6 = 1 nm, type II: b — 0.5 nm). The lowest exciton levels are shown in 
a) and d) (for details see text). The spin-dependent Zeeman energy is not included. In b) and c) the corresponding oscillator 
strengths are plotted. The lowest panel c) and d) show the second derivative of the ground state energy which quantifies 
the X-ABE oscillations. 



The period can be estimated from the second line of 
Eq. ([T^ . assuming that due to the confinement both (re) 
and (r/i) are almost constant, as 

2^1 1 

which is Bp — 20.6 T in the present case. The minima 
of the oscillatory component of the energy are found at 
B = jBp, j being integer, and maxima at _B = jBp/2. 
Consequently, the minimum of the second derivative is 
found at around half of the oscillation period in all cases 
{Bp/2 = 10.3 T). 

The oscillation amplitude is reduced for the eccentric 
nanoring. Here, the non-uniform ring width tends to 
push the exciton wave function into the broader part, 
thus weakening its ring topology which is a necessary 
ingredient for X-ABE oscillations. Note that for b > r2 — 
ri the confinement potential reduces to a banana-shaped 



quantum dot, and the ring topology is lost completely. 



B. Type II nanoring 

Secondly, we investigate a type II nanoring consisting 
of InP in the ring and Gao.51Ino.49P outside, embedded 
into an AlAs barrier along z. Strain effects are taken into 
account only in so far as the hole is always found inside, 
r/i < ri . A full inclusion of strain would modify the exact 
potential profile but not the 0(2) or Di symmetry. 

For the concetric type II ring. Fig. [5ji illustrates the 
crossings among states with different quantum number 
La (black solid La = 0, dashed La = ±1, and dotted 
La = ±2). This crossing resembles more the ABE for 
individual carriers and can be traced back to the much 
reduced exciton effect in the present type II nanoring. In 
order not to overload Fig. Oi; some higher states have 



been omitted. If the symmetry is reduced to Di (red 
[grey] curves), all states are mixed similar to the type 
I nanoring. Consequently, all crossings become anti- 
crossings and the energy dependence on the magnetic 
field differs not much fi'om the type I nanoring, with re- 
duced X-ABE oscillations (Fig. [Sf). 

The difference between type I and II nanorings be- 
comes apparent when the oscillator strength is studied as 
depicted in Fig. [Sj3. The oscillator strength correspond- 
ing to different exciton states changes with the magnetic 
field as the character of the exciton wave function changes 
itself: The main component of the ground state wave 
function isL = Oati? = OT and shifts to L = 1 at 
B = 10 T. This results in a decrease of the oscillator 
strength which is transferred to the first excited state 
(and increased due to the larger electron- hole overlap). 
Later on it is transferred to the second and higher excited 
states where different quantum numbers L mix strongly. 
Even though there is not any strict rule for the oscil- 
lator strength conservation. Fig. [3^ suggests that it is 
approximately valid for the lowest exciton states. The 
rapid decay of the ground state oscillator strength with 
increasing magnetic field has interesting implications for 
the exciton kinetics (Section IV). 

Unlike the type I nanoring, the oscillation period 
changes with b (Fig. [Sf). This is also accompanied by 
a decrease of the oscillation amplitude due to the loss 
of the ring topology with increasing b. The oscillation 
period is now determined by the periodic change of the 
ground state main component from L to L -I- 1, thus by 
the center-of-mass motion. The exciton relative motion 
with the period Bp in Eq. (PT|) plays only a minor role 
for every value of 6 ^ 0. The value of the center-of-mass 
period can be estimated from the last term of Eq. ^ 
assuming again strong radial confinement. 
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In the present case, taking i?g = {r^) = 12.3 nm and 
Rh = {fh) = 4.4 nm we obtain Bp,2 ~ 9.9 T, which agrees 
well with the observed period in Fig. [5f. 

In Fig. [S] correlated one-particle densities arc plotted. 
They are defined as conditional probability to find a par- 
ticle, either electron or hole, in the exciton while fixing 
the coordinates of the other particle at a certain position 
marked by a large dot. For the formal definition, see 
Eqs. (|A11 IA2[) . This concept has been used also for an- 
alyzing the two-electron Wigncr molecule in a quantum 
dot.— An inspection of Fig. [6] reveals how the loss of the 
circular symmetry in the eccentric nanoring modifies the 
wave function. By fixing the hole position at the right 
(wide) side of the ring, most of the electron density is 
found there as well (Fig.[6li), which resembles the situa- 
tion in a concentric ring (not shown). The same is true 
for the hole (Fig. [6];). On the other hand, by fixing either 
the electron or the hole on the opposite narrow side of 
the ring, the picture changes. Due to the strength of the 
ring confinement and the weakness of the Coulomb in- 
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FIG. 6: Correlated electron and hole densities of the exciton 
ground state at zero magnetic field for a non-circular type II 
nanoring (InP/GalnP, ri = 8nm, r2 = 16 nm, b = 2nm). 
The black dot indicates the chosen coordinates of the other 
particle. 



teraction in type II structures, the electron is still found 
on the right where the confinement energy is minimal 
(Fig. ^p) . Since the electron density is very small on the 
opposite side, the Coulomb correlation seen by the hole 
is tiny, and an almost circular - i.e. one-particle like - 
hole density is found in Fig. [6]i. 

Finally, let us briefly consider the case of a free 
electron-hole pair which has been investigated for type 

I nanorings in Ref. [13. Formally, we switch off the 
Coulomb interaction Vc in the Hamiltonian Eq. ([T]) . The 
free electron-hole transitions shown as dashed curves in 
Fig. [7^ are dominated by the electron level since the hole 
part refers to an almost fixed angular quantum number 
due to the small effective hole radius of 4 nm. Therefore, 
the first levels can be characterized by the electron quan- 
tum numbers = 0, — Their degeneracy at zero 
magnetic field is lifted when switching on the Coulomb in- 
teraction (solid curves), similarly to Fig. [SJl. The overall 
down shift in energy is a measure of the exciton binding 
energy. 

The second derivative displayed in Fig. [7)3 shows that 
the diamagnctic shift (being proportional to the relative 
distance) of the free electron-hole pair is larger than in 
the exciton. The reason is the Coulomb attraction in the 
exciton which brings electron and hole closer together. 
The oscillation amplitude, however, is getting stronger 
for the exciton. We conclude that the ring topology of the 
electron part is stabilized due to the Coulomb attraction 
with the hole. Therefore, surprisingly in the specific type 

II example, the ring eccentricity is felt not as strong for 
the exciton compared with the free e-h- pair. 
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FIG. 7: (Color online) Calculated exciton (solid) and free electron-hole pair (dashed) properties in non-circular nanorings of 
type II (InP/GalnP) with ri = 8nm, r2 = 16 nm, and 6 = 2 nm as a function of magnetic field. The ground state in black, the 
first and the second excited state in red (dark grey) and green (light grey) are shown in a) (the spin-dependent Zeeman energy 
is not included). In b) the second derivative of the ground state energy is depicted. 



To summarize this section, we have demonstrated that 
the amplitude of X-ABE weakens when going from cir- 
cular to non-circular symmetry since the wave function 
tends to lose its ring topology. However, if the asymme- 
try is weak and electron and hole are spatially separated 
as in type II structures, the ground state still exhibits 
oscillations in its energy, while being always bright and 
observable by optical means. 



IV. KINETIC EQUATIONS 

For calculating the photoluminescence emitted from 
the nanoring, we need to know the occupation Na of 
each exciton states. For the linear density regime and in- 
coherent excitation, the relevant set of kinetic equations 
has been derived in Ref. 



~dt 



= 9o 



+ ^7,3o)A^o.(23) 



Here, ga is a state dependent generation term which 
stands for the last term in a chain of optical phonon 
emission events after optical intcrband excitationfS^ The 
radiative decay rate Tq, of a localized exciton state con- 
tains emission into both TE and TM polarization and is 
found proportional to the squared oscillator strength)^ 



3 h^ci 



/a I 



(24) 



where nji is the refractive index and Eg the band gap. In 
Eq. da is a phenomenological non-radiative decay 

rate, representing processes as e. g. exciton annihilation 
via impurities, escape into the wetting layer, or Auger 
processes. These processes are not treated explicitly in 
the present theory. 

The acoustic phonon scattering rates 7q,^ are defined 



7a/3 



27r 
q 



nB{'hujq)5{Ep - Ea + TlOJq)] 



(25) 



where nB{hLOq) is the Bosc-Einstcin distribution of acous- 
tic phonons with dispersion fiuiq = hsq {s - sound veloc- 
ity). They obey the relation of detailed balance between 
in- and out scattering of a given state. 



(_B„--Bfi)/fcf3T 



(26) 



with the phonon (i.e. lattice) temperature T. Strictly 
speaking, lattice vibrations in nanostructures differ from 
the respective bulk ones;^ but for the present purpose 
this refinement is of minor importance4^ The evaluation 
of the exciton-phonon matrix elements for deforma- 
tion potential scattering is given in the Appendix. 

In the exciton scattering with acoustic phonons we 
have taken into account only the spin diagonal part. In- 
deed, there is non-diagonal scattering (" spin-flip" ) as well 
due to the lack of the inversion symmetry in nanostruc- 
tures - Rashba efFect4i Under linearly polarized excita- 
tion, the initial generation is spin independent. As- 
suming that in the final stage of thermalization spin-flip 
processes can be neglected, we expect an equal occupa- 
tion of spin up and down bright states. This is supported 
by recent experiments showing almost no difference in in- 
tensities between cr"*" and a~ polarized lines4^ 

The kinetic equations are solved numerically for the 
steady state situation dNa/dt = 0. The resulting occu- 
pations Na enter the photoluminescence spectrum I{E), 



I{E) = ^r„iV, (?(£;-£;„). 



(27) 



Note that the linear absorption spectrum D{E) does not 
depend on the occupations. 



D{E) = ^r„<5(i; -£;„). 



(28) 



9 



Circular Non-circular 







1 1 1 1 1 1 





J I I I I I I U I I I I l_ 

5 10 15 20 25 5 10 15 20 25 

B (T) B (T) 



FIG. 8: Calculated absorption spectra of a circular (6 = nm) 
and non-circular {b = 0.5 nm) type II nanoring (InP/GalnP, 
ri = 8 nm, r2 = 16 nm) including the spin contribution 
Eq. ([2]). All spectra are Gauss broadened with variance 
CT = 0.1 meV and displayed using a linear gray scale. 

Again, a constant prefactor has been omitted. 

Summing Eq. (|23p over all states a, the phonon scat- 
tering terms cancel, and the following conservation law 
for the total pump rate P is found 

P = ^5„ = ^(r„-Hda)iVa. (29) 

a a 

We will exploit this relation in Section V in the discussion 
of PL quenching, i.e. a decrease of the PL intensity with 
S-field. 



V. ABSORPTION AND 
PHOTOLUMINESCENCE 

Since in type I nanorings, the exciton energies and 
oscillator strengths do not depend much on the mag- 
netic field, interesting effects like PL quenching are not 
to be expected. Therefore, we focus here on a type II 
system and investigate InP/GalnP nanorings with radii 
ri = 8 nm, r2 = 16 nm. Results for an eccentric nanoring 
(c) are compared with a sligthly eccentric one (n), hav- 
ing a center displacement 6 of 0.5 nm as used in Fig. [51 
right panel. The radiative rates of the lowest state a = 
arc rather small (rg = 0.042 ns""'^ and Tq = 0.027 ns~^ at 
S = T) due to the tiny overlap between electron and 
hole in the wave function, as specific for any type II struc- 
ture. Similar small radiative rates have been calculated 
for spatially indirect excitons in coupled quantum wells 4^ 
The exciton-phonon scattering rates from the first ex- 
cited state down to the ground state are ^'^■y = 102 ns~^ 
and 7q\ ~ 138 ns~^ at i? = T. These numbers clearly in- 
dicate that the exciton-phonon scattering dominates the 
kinetics. 

The absorption spectra of the circular and non-circular 
nanorings plotted in Fig. [8] are calculated for linearly po- 
larized light where both spin components Eq. ^ with 
M ~ ±1 are present. Therefore, all lines appear as Zee- 
man splitted doublets. 
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FIG. 9: Calculated photoluminescence spectra of a circular 
(a and c) and a non-circular (b and d) type II nanoring as 
in Fig. [8l Exciton occupations are calculated from the ki- 
netic equation Eq. (|23|l . fixing the phonon temperatures as 
indicated. Without non-radiative decay, da — 0. 



The absorption spectra plotted in Fig. [5] show pro- 
nounced differences between both nanorings: In the con- 
cetric nanoring only two doublets having quantum num- 
ber La = are visible (Fig. [Hi). Above B = 5T, the 
lower state is not any longer the ground state (compare 
Fig. [5ji). In the eccentric ring the oscillator strength is 
transferred from the ground state to higher states with 
increasing magnetic field (Fig. [8)d) as already discussed 
in detail in Section III-B. Although the oscillation of the 
lowest bright state can be hardly seen in both cases, their 
second derivatives reveals them clearly (see Fig.^). The 
excited states exhibit much stronger oscillations since 
here the Coulomb attraction between electron and hole 
acts much less. Both periods agree well with Eqs. ((2T|) 
and ([22]) . respectively. 

In view of an easier interpretation of the calculated 
spectra, we envisage in Figs. |9] and [10] a PL detection 
with circularly polarized light (T+. Thus, only the lower 
Zeeman splitted lines are seen. For solving the kinetic 
equations, nine exciton states were taken into account. 
First, we investigate the case without non-radiative de- 
cay, da — (Fig. [5]). Due to the large energy differ- 
ence between the first and the second bright state, only 
the lowest one is visible for both temperatures (Fig. [9^ 
and c). The results for the eccentric ring are not so 
simple. Although the oscillator strength of the ground 
state goes down appreciably (at B = 25 T it is less than 
10^^ of its value at B = OT), the luminescence line at 
T = 4K (Fig. [9)3) has almost constant intensity, and 
no line quenching is seen. This is a consequence of the 
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FIG. 10: Calculated photoluminescence spectra as in Fig. [O] 
but here with a state independent non-radiative rate of da = 
10ns"\ 



conservation law Eq. ([29|) and will be discussed below 
in more detail. At elevated temperature, other lines are 
seen as well, but their intensities have a different B de- 
pendence compared to the absorption, which signals the 
role of exciton occupation (Fig.[5Ji). 

Second, non-radiative decay is included in Fig llOl For 
the discussion, three rates respectivelly times are impor- 
tant: (i) the radiative rate (tens of nanoseconds), 
(ii) the exciton-phonon scattering rate 7q,/3 (several pi- 
coseconds), and (iii) the non-radiative rate da which is 
a phcnomcnological input here. With the assumption 
= 1/(100 ps) wc use a value which dominates over 
the extremely small radiative rate, but is well below the 
phonon scattering rate. Wc start again with the discus- 
sion of the circular ring where the situation is simple. 
Due to the change of the ground state from bright to 
dark one at i? > 5 T, these excitons decay predominantly 
non-radiatively, which results in a distinct line quenching 
(Fig. [TUk and c). On the other hand, for the non-circular 
nanoring more lines are seen in the PL, in particular at 
low temperatures (Figs.lTOb and d). Obviously, complete 
equilibration with the lattice temperature is no longer 
reached. For the unrealistic case of non-radiative rates 
being stronger than the phonon scattering, the exciton 
occupation is simply given by the ratio between pump 
rate ga and decay da- Keeping both quantities constant, 
the exciton occupation gets constant, and the PL spec- 
tra would coincide with the absorption spectra. Although 
this extreme limit is not reached with the actual param- 
eters, the tendency is clearly observable in Fig. [TOl 

Let us now concentrate on the integrated photolumi- 



nescence / which follows from Eq. ((27|) as 



= > TaNa. 



1 = J dtoliu) = 



(30) 



Its relation to the pump rate P will be called photolumi- 
nescence yield Y , 



(31) 



Without non-radiative decay {da = 0), the conservation 
law Eq. p9|) gives immediately Y = 1 independent of B 
- each excited exciton (or in general electron-hole pair) 
decays into one emitted photon (full curve in Fig. [TTk 
and b). 

Things change a lot if non-radiative decay channels are 
included, da > 0. Due to the extremely small radiative 
rates, these processes can even dominate the exciton de- 
cay, leading to a yield much below unity fFig.fTTh and d). 
Since the ground state of the concentric ring is getting 
dark at around B = 5T, a steep decay of the yield fol- 
lows. In the eccentric nanoring the yield goes down not 
as abruptly since the oscillator strength of the ground 
state decays more slowly. The slight oscillations seen in 
the full curve of Fig.fTTH are related to the changing level 
distances which influence the individual phonon scatter- 
ing rates. We conclude that the quenching of the total 
PL is intimately related to non-radiative processes. The 
quenching is not as dramatic at elevated temperatures 
since higher exciton states contribute more and more to 
the total emission (not shown). 

In all cases studied here, the phonon scattering rates 
are dominant. Therefore, it can be expected that the 
occupations of the different exciton states deviate not too 
much from equilibrium which in the present low-density 
case, in accordance with Eq. ([26]), is characterized by the 
Maxwell-Boltzmann (MB) distribution 

Na = Cexp{-PEa) (32) 

Within this approximation, the PL 



with P = l/ksT 
yield is given by 



Y. 



MB 



-/3-Eq 



(33) 



and shown as dotted curves in Fig. IIH indeed not very 
much different from the full calculation. This is the right 
place to discuss Ref. [l^ which was the first attempt to 
calculate the exciton PL for a nanoring. Concerning 
the exciton kinetics, they (i) have assumed a Maxwell- 
Boltzmann distribution for the excitons and (ii) nor- 
malized the integrated PL to constant exciton density 
iV = Na - and not to pump rate P - resulting in 



Y. 



MBN 



-/9-E„ 



y 



(34) 



This normalization makes a pronounced difference when 
non-radiative processes arc absent (dashed curve in 
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FIG. 11: (Color online) The photoluminescence yield Y for the circular (a and c) and non-circular (b and d) nanorings of Fig.|8] 
at T = 4K and two values of the non-radiative decay rate da — a) and da = 10 ns"^ b). The full solution Eq. (|31|l (Full, 
solid) is compared to the assumption of a Maxwell-Boltzmann distribution with constant pump rate Eq. ([32} (MB, dotted) 
and with constant exciton density Eq. (|34p (MEN, dashed) is plotted. 



Fig. [TTk). but gives nearly identical results with Eq. (1551) 
when these dominate (dashed curve in Fig. [Tib). Fmbn 
has been normalized to Ymb at _B = T. 

While this discussion refers to PL quenching, in Ref. [T3 
PL blinking was proposed, too: The integrated PL goes 
down and up in dependence on B, since the ground state 
switches between bright and almost dark behavior. The 
authors have used a ring model with zero width and 
rather close radii for electron an hole. The latter seems to 
be decisive for the blinking effect to occur. For the type II 
nanoring studied in the present work, such a blinking can- 
not be expected since the average electron and hole radii 
are rather different {Re = 12.1 nm vs. Rh = 4.4 nm). 
Moreover, the blinking effect predicted in Ref. [13 has 
been questioned recently by noting that a better account 
of the Coulomb interaction is neededfi^ 



VI. CONCLUSIONS AND OUTLOOK 

We have derived the Hamiltonian for excitons in nanor- 
ings with finite width and arbitrary symmetry. Sorting its 
matrix elements according to the symmetry point group 
of the confining potential we have shown how the dif- 
ferent wave function components are coupled. Two pro- 
totype systems of type I (GaAs/AlGaAs) and type II 
(InP/GalnP) have been considered. The numerical in- 
vestigations showed that: (i) The oscillation amplitude of 
the exciton Aharonov-Bohm effect decreases when going 
from circular to non-circular symmetry due to the addi- 



tional localization of the exciton. Moreover, the oscilla- 
tion period in type II nanorings changes from relative- 
motion induced to being determined by periodic changes 
of the center-of-mass wave function, (ii) The exciton 
ground state in the non-circular (eccentric) nanoring re- 
mains always optically active since the total angular mo- 
mentum is no longer a good quantum number. However, 
its oscillator strength can be extremely small in type II 
nanorings. 

Further, we have investigated the exciton kinetics 
within a model which includes acoustic phonon scatter- 
ing and radiative and non-radiative decay. Our study 
of a slightly eccentric type II nanoring has revealed that 
the oscillations of the excited states are clearly visible 
and that the amplitude of the ground state oscillations is 
weaker than in the concentric ring. Calculating temper- 
ature dependent occupations of exciton states we could 
show how the appearance and disappearance of single 
lines in PL are related to their counterparts in absorption. 
The presence of non-radiative decay is decisive for the 
quenching of the integrated photoluminescence. If the 
exciton-phonon scattering dominates over all decay rates, 
a simplified description using a Maxwell-Boltzmann ex- 
citon distribution works reasonably well. 

One unexpected finding is that a small violation of 
the circular symmetry improves the observability of the 
X-ABE since here all states are optically active, and anti- 
crossings periodic with B can be easily seen in the second 
derivative of the ground state energy. Non-radiative de- 
cay channels allow to sec more lines in PL, thus improv- 
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ing the observability of X-ABE for excited states even 
at low temperatures. Therefore, a slightly asymmetric 
nanoring of type II seems to be the best candidate for 
the experimental confirmation of the exciton Aharonov- 
Bohm effect in linear optics, i.e. absorption or photolu- 
minescence. From a practical point of view, an exactly 
circular nanoring would be rather the exception than the 
rule, given the uncertainties of nanostructurc growth on 
patterned substrate, not to think of self-organized ring 
formation. 

Quite recently, the possibility of implementing a flux 
qubit in small non-superconducting rings has been dis- 
cussedi^ Since a persistent current due to excitons can 
be initiated and controlled optically^^ it can be specu- 
lated that an exciton qubit in a nanoring can be formed 
whose function rests upon the exciton Aharonov-Bohm 
effect. 
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APPENDIX A: TECHNICAL DETAILS 



fective parameter which has been fitted to give reasonable 
agreement with the quantum well Coulomb potential. 
The phonon matrix elements are defined as^ 

(^Dcexp(iqre) - cxp(iqr,,)^ $/3(re, r,i), 

(A3) 

where pM is the mass density, V the sample volume, and 
Dc {Dy) the deformation potential for electron (hole). 
Using the single sublevel approximation and the expan- 
sion of the wave function Eq. ([5]), we obtain in cylindrical 
coordinates for the three-dimensional momentum q 

(A4) 

introducing the state dependent overlap functions 
^L^L'^'^\l\\) ^^'^ ^ z-dependent contribution 



The correlated one-particle densities are defined as 
conditional probabilities 



^'a r,rh,<l>^ (fi - <l>h, $ 



<Pe + V 



<l>h + (P 



(Al) 

2 



(A2) 



where the coordinates of the other particle are fixed at 
(r. If) . Note that this definition differs from the previous 
one in Ref. [l^ where an additional integration over r was 
performed. 

The averaged Coulomb potential is for convenience ap- 
proximated by^ 



Voire, rh,( 



1 ■ u fbc'L, 

-arcsmh 



where r = ■\/r^ + rf^ — ^r^rh cos(0), and is the width 
of the quantum well (4 nm in our case), be = 0.6 is an ef- 



/•oo 

I 

><'"zzp.kz^,L',/3(^e,?'h) JL-L'{q\\re(h)); (A5) 

Ka{q.) = [ dzvl{z)e-'^'' . (A6) 



Jl-L'{x) arc Besscl functions of the first kind, Va{z) is 
the confinement wave function in z-direction, and the 

symmetry Sl'^l^'''\q\\) = slT^'^\q\\) holds. In order 

+1 



to get the scattering rate 7^/3, the matrix element 



squared has to be integrated over q. The integration 
over q^ is performed using the energy conserving delta 
function in Eg. [25l which gives |Aq,/3| = Tisq with the en- 
ergy difference between states A^^ = — Ep . Defining 
the overlap sum 



W^'piqii) 



LAIN 



the final expression is obtained as 
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lap 



27rh^s^t 



dq\\ 911 



PM Jo - qpq^ 



(A8) 



where qz 



Note that both phonon emis- 



sion and absorption processes (first and second term in 
Eq. (P5)l . resp.) are included here, since A^^ can have 
both signs. 

In order to simplify further, an approximation intro- 
duced in Ref. IH is adopted here, too. The strongest 
confinement is found in the growth direction z which al- 
lows to put K'^{qz) ^ K^iO) ~ 1- Further, a rapid decay 
of W^^(g||) is assumed, well before the integration limit 
q is reached. Then, the integral in Eq. (jA8[) can be ap- 
proximated by 



dq\\ q\\ 



1 - qi/q' 



dq\\ q\\ 



which allows to integrate over g|| analytically usingsi 
/ dq\\ g|[ JN{raq\\) JN{rbq\\) = —S{ra - n) . 



The general expression Eq. (|A8p reduces to 

7.. = "^^fff"^ X (A9) 

introducing the following abbreviations: 

/>oo 

I Jo 2 ' ' 

^LL^ir) = XI dr' r' ui^lA^' ,r)u^^L^^^,^^{r' ,r) . 



The detailed balance Eq. (|26p can be checked in Eq. (jA9 
quite directly. 
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